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We show that guiding of optical signals in chains of metal nanoparticles is subject to a surprisingly 
complex dispersion relation. Retardation causes the dispersion relation to split in two anticrossing 
branches, as common for polaritons. While huge radiation losses occur above the light line, just 
below the light line the micron-sized loss lengths are much longer than expected. The anticrossing 
allows to create highly localized energy distributions in finite arrays that can be tuned via the 
illumination wavelength. Our results apply to all linear chains of coupled resonant scatterers. 

PACS numbers: 42.25.Fx,42.79.Gn, 78.67.Bf, 71.45.Gm,73.22.Lp 



Ordered arrays of optically driven metal nanospheres 
may be used to transport optical signals in structures 
that are much smaller than the wavelength of light 
Q IE Q 1^ ' Coherent transport occurs for frequencies 
close to the plasmon resonance of a single nanosphere, 
and is mediated by dipole-dipole coupling between reso- 
nant particles on a chain. The transverse width of the 
guiding structure is only limited by the particle diame- 
ter, which is an order of magnitude smaller than the opti- 
cal wavelength. Metal nanoparticle arrays may therefore 
allow to overcome current limitations in the miniatur- 
ization of integrated optical devices. Although several 
experimental and theoretical studies [E l^b i^i i2n^] 
confirm that coupling effects occur in metal nanosphere 
arrays, no complete picture for the dispersion relation 
and propagation loss has emerged. 

Recent studies indicate that radiation damping, re- 
tardation effects and long-range coupling, dramatically 
affect the loss and dispersion 0, |g- This necessi- 
tates a complete revision of the original theoretical stud- 
ies 0, 0, 0, . In this paper we calculate the dispersion 
and loss including all these effects for finite and infinite 
arrays. Remarkably, retardation causes the dispersion to 
split into two anticrossing branches for transverse modes. 
Modes with wave vectors slightly larger than those in the 
embedding medium have the lowest loss, and allow trans- 
port over micrometer distances, well in excess of previ- 
ous estimates 0. We further show that the anticrossing 
gives rise to a strongly localized optical response in finite 
arrays that is very sensitive to the incident wavelength 
around the anticrossing range. As our calculations are 
not specific to metallic nanoparticles, our results apply 
in general to periodic arrays of coupled resonant dipoles, 
in nano-optics as well as in, e.g., atomic arrays interact- 
ing with radiation j^, llfl . 

We consider finite and infinite arrays of equally spaced 
metal nanospheres of radius a, spaced by a center-to- 
center distance d. A coupled point-dipole approximation 
is well suited to describe the electromagnetic response of 
such chains 0, S ^| ■ We use the electric field gener- 



ated by a single dipole pe , oscillating with frequency 
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where f is the unit vector pointing from the dipole to 
the field point at distance r, e = n^eo is the permittiv- 
ity of the homogeneous medium in which the dipole is 
embedded, and v = c/n is the corresponding speed of 
light. This form of the electric field fully takes into ac- 
count retardation effects. The induced dipole moment 
on a nanosphere equals its polarizability a(uj) times the 
local electric field, which is composed of the driving field 
g(cxt) g^^^ ^Yie fields of all the other nanospheres. For the 
n*'' nanosphere in a linear chain pointing along f we find 
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For N dipoles, Eq. represents N inhomogeneous cou- 
pled linear equations for the dipole moments p„. 

The resonant material response of the metal 
nanospheres is captured in the frequency dependence of 
the polarizability a{uj). Within the Drude model, the 
polarizability of a nanosphere is Lorentzian 



aDrude(w) = 47re- 
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Throughout this work, we use values 7 = 8.25 • 10^'' s~-^, 
and Losp = 4.64-10^^ rad/s appropriate for silver particles 
in glass {n = 1-5}, as derived from tabulated optical con- 
stants in Ref. [l3|- To maintain the energy balance be- 
tween extinction, scattering and absorption it is essential 
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FIG. 1: Excitation by a plane wave along the array causes 
a localized optical response in finite plasmon arrays (10 Ag 
particles in glass, radius 30 nm, d — 75 nm). Depending on 
driving frequency, the Ohmic dissipated power (at 1 V/m in- 
cident field) is either concentrated on the first (solid curve) , or 
last sphere (dashed) encountered by the incident beam. Ver- 
tical dashes indicate the single particle resonance frequency. 



to include radiation damping in the polarizability of each 
nanosphere by setting the polarizability to j7|, |8,, .12.] 
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The importance of retardation effects in the optical 
behavior of metal nanoparticle chains is evident from 
Figure ^a), where we consider the optical response of 
a chain of 10 nanoparticles of radius 30 nm and spacing 
75 nm illuminated by a plane wave incident along the 
array. We plot the Ohmic dissipated power (oc |PnP) 
for the first and last sphere in the array as a function of 
the driving optical frequency. In agreement with results 
reported by Hernandez et al. [l^ . we find a large asym- 
metry in the response of the array, which can be tuned 
via the driving optical frequency. For frequencies below 
4.3 • 10^^ rad/s the backmost nanosphere is preferentially 
excited (dashed line in Fig.^, i.e., the sphere that is en- 
countered last by the excitation beam. In contrast, this 
sphere is hardly excited for higher frequencies, for which 
the frontmost sphere is strongly excited (solid curve). 

Naively, one expects such strong coupling between the 
incident wave and the plasmon chain to occur when the 
plasmon dispersion relation intersects that of the em- 
bedding medium ('light line'), i.e., when the incident 
field is phase matched to the mode in the array. For 
metal nanoparticle chains a quasistatic model is often 
used 1^ Q , that is found from Eq. ((SJ for an infinite 
number of dipoles, in the limit 7 = 0, c ^ oo and in 
absence of an external driving field. For the transverse 
modes that the incident plane wave considered in Fig- 
urencan couple to (pm -f = 0), the quasistatic dispersion 
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where k is the wave vector. This dispersion relation is 
plotted in Figure A) (red solid line). Evidently, the 
light line (dashed) intersects the quasistatic dispersion 
relation at a frequency of 4.6T0^^ rad/s that does not 
correspond at all to the frequencies around 4.3-10^^ rad/s 
for which the complex response in Figure ^ occurs. Fur- 
thermore, our exact calculations show that the asymmet- 
ric response in Fig. ^ of arrays under plane wave illumi- 
nation shifts further to the red away from the plasmon 
resonance as the spacing between particles is increased 
(data not shown). This shift is in sharp contrast to the 
rapid decrease of the bandwidth of the quasistatic 
dispersion relation with increasing particle spacing, that 
is caused by the strongly reduced overlap of particles with 
the fields of their neighbors. We conclude that the asym- 
metric response in Fig. ^ points at the strong effect of 
retardation on the array dispersion. Retardation affects 
the phases of individual nanospheres, and thereby the 
constructive or destructive interference along the chain. 

Damping and radiation losses preclude the existence 
of a dispersion relation which gives solutions of the form 
Pm oc e""'^'' for real frequencies. However, complex fre- 
quency solutions exist, corresponding to the dispersion 
of damped modes. Recently, several authors have at- 
tempted to determine the complex dispersion relation of 
plasmon chains. Weber and Ford 7] calculated the com- 
plex eigenfrequencies for finite chains of N dipoles, which 
are frequencies for which the matrix coupling the dipoles 
in Eq. (01 is singular. A unique wave vector can be as- 
signed to these N modes based on their mode profile ■ 
The resulting set of solutions is expected to be an accu- 
rate but discrete approximation to the infinite chain dis- 
persion relation, that includes effects of retardation, and 
ohmic and radiation damping. Ref. and open symbols 
in Fig. IHA) show that the resulting dispersion, defined 
via the real part of the complex frequencies, is dramati- 
cally different from the quasistatic result, especially near 
the light line. Unfortunately, the discrete sampling that 
is inherent in this approach does not allow to distinguish 
between the two distinct scenarios, in which the disperion 
relation either has a polariton form with two branches 
that anticross at the light line, or is a single continuous 
dispersion relation. 

To resolve this ambiguity, we calculate the infinite 
chain dispersion by inserting p,„ oc e*™'^'' in Eq. and 
setting the driving field to zero. Citrin Q recently re- 
ported that the sums, which exhibit poor convergence for 
frequencies with negative imaginary parts, can be evalu- 
ated in terms of poly logarithms |l5| Li„, resulting in an 
implicit dispersion relation for transverse modes of the 
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with 
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In Ref. j^] this solution was used to approximate the infi- 
nite chain dispersion perturbatively by evaluating E(w, fc) 
at the resonance frequency ujsp- The result, plotted as 
the dotted line in Fig. |2IA), shows that this approxi- 
mation differs significantly from the quasistatic result 
and has a deep minimum where the light line crosses 
the quasistatic dispersion. This minimum is due to 
the 1/r term in the dipole field, which contributes the 
term [Lii(e*("/''-'=)'^) -h Lii(e*('^/^+'=)'^)] with a logarith- 
mic singularity at the light line. As u; is far from losp 
the validity of this perturbative approach is limited. 
Therefore we have numerically solved the full dispersion 
relation to find the complex frequencies uj corresponding 
to real wave vectors. 

FigureEfB) (solid red dots) shows the dispersion of the 
real part of the frequency for an infinite chain of a = 30 
nm radius silver particles spaced by d = 75 nm, embed- 
ded in glass {n = 1.5). First of all, our result confirms 
that the deviations of the discrete finite chain result in 
Fig. ElA) from the quasistatic result for infinite chains 
are not due to chain length effects, but to retardation, 
radiation damping and ohmic damping. More surpris- 
ingly, we find that the dispersion relation for transverse 
excitations has two branches that anticross at the light 
line, indicating polariton behavior of the coupled photon- 
plasmon chain system. Regarding the upper dispersion 
branch, we note that this branch does extend upwards 
along the light line, though our numerical algorithm does 
not identify any clear roots as loss rates become compa- 
rable to and exceed the optical frequency. Indeed, the 
dispersion branch above the light line is accompanied by 
huge damping, as gauged by the imaginary part of the 
frequency plotted in Fig. EJC). Damping times for this 
branch are shorter than about 1 optical cycle, and radia- 
tive losses diverge as the light line is approached. 

The existence of two dispersion branches that display 
an avoided crossing is consistent with, but could not be 
surmised from the finite chain result obtained by Weber 
and Ford and is not represented by the perturbative 
treatment by Citrin [i| (dotted line in Fig. I^A)). That 
two disconnected branches appear may seem surprising 
in view of textbook treatments of, e.g., the dispersion 
of surface plasmon polaritons. Damping causes separate 



branches to merge when real frequencies but complex k- 
vectors are considered 0|. The dispersion of complex 
frequencies for solutions with real wave vectors, how- 
ever, displays an avoided crossing, also in the presence 
of damping. A treatment in terms of complex frequen- 
cies, i.e. normal modes that decay in time, is the natural 
extension of the normal mode analysis for finite arrays. 

We now turn to the relevance of the exact disper- 
sion relation derived here for proposed applications of 
metal nanoparticle chains as ultrasmall optical waveg- 
uides. Both the group velocities and the damping times 
in the present calculation are profoundly different from 
those found in the quasistatic approximation ^IBj- 
In Figure|2IC, inset) we consider the propagation length, 
defined as the product of the damping time in Fig. |2IC) 
and the group velocity that we derive from Fig. H^B). 
Based on the quasistatic model, the longest propaga- 
tion lengths are expected near the center of the Brillouin 
zone [2I, y. Such predictions as well as extrapolations 
from k = measurements need to be completely re- 
vised due to the polariton splitting found here. For trans- 
verse excitations. Fig. I3C) shows that modes above the 
light line are strongly damped, with damping times com- 
parable to the optical period, and propagation lengths 
less than the interparticle spacing (solid red dots). For 
modes just below the light line {k > 0.37r/ci), however, no 
radiative losses occur. Damping rates of around 10^'^ s~^ 
and group velocities around 0.3c provide a frequency win- 
dow for which propagation lengths above 5 fim are pos- 
sible. 

Applications that require guiding can benefit from lon- 
gitudinal modes, for which the induced dipoles point 
along the array. As was also found for finite arrays by 
Weber and Ford, retardation effects cause the exact dis- 
persion relation (open squares in Fig. |21[B)) to have a 
doubled bandwidth compared to the quasistatic result 
(dotted curve). As longitudinal modes do not couple to 
plane waves that propagate along the array, the intersec- 
tion with the light line is not associated with diverging 
loss (Fig.[2IC), open squares) and no splitting of the dis- 
persion relation occurs. In general, guiding benefits from 
the fact that the group velocity is large over a much wider 
bandwidth. Below the light line a large frequency window 
occurs for which the propagation length exceeds 1 /xm. 

In conclusion, we have presented the dispersion rela- 
tion of infinite metal nanoparticle chains, fully taking into 
account the effects of Ohmic damping, radiation damping 
and retardation. For transverse modes strong coupling to 
plane waves propagating along the array causes an un- 
expected splitting of the dispersion into two branches, as 
common for polaritons. Propagation lengths just above 
1 /im for modes below the light line may allow plasmon 
chains to be used as very small guiding structures for 
nanoscale energy transport. With respect to the asym- 
metric optical response in Figure^ it appears that the lo- 
cal excitation intensity set up by an incident plane wave. 
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as gauged by the Ohmic dissipation per sphere, can be 
used as a a sensitive fingerprint of the plasmon chain po- 
lariton anticrossing. The sharply defined frequency at 
which the front-backward asymmetry changes direction, 
and the peak absorption frequencies for the front- and 
back sphere provide a measure for the location and mag- 
nitude of the splitting in the dispersion relation. Such 
steep frequency edges are sensitive to the refractive in- 
dex of the medium surrounding the array, and can be 
useful for, e.g., nanoscale sensing applications. Further- 
more, the arrays allow to create strongly localized energy 
distributions, that can be varied by tuning the incident 
plane wave around the anticrossing in the dispersion rela- 
tion. Such tunable localized energy distributions are use- 
ful for, e.g., locally enhancing nonlinear interactions, and 
nanoscale photolithography with visible wavelengths. Fi- 
nally we note that these phenomena should occur for any 
chain of coupled resonant dipoles, which also includes 
systems like quantum dot arrays, or atomic systems 
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Exact results 




k (units n/d) 



FIG. 2: (Color online) (A) The quasistatic dispersion (solid 
red curve) for transverse modes intersects the light line 
(dashed) slightly below the single particle resonance (hori- 
zontal dashes). Blue open circles show the 10 normal mode 
frequencies for an array of 10 particles, according to Ref. 0, 
which agree neither with the quasistatic dispersion, nor with 
the perturbative result of Ref. @] (dotted curve). (B) Red 
solid circles (green open squares): real part of the frequency 
for roots of the exact infinite chain dispersion for transverse 
(longitudinal) modes. For transverse modes two branches ap- 
pear that anticross at the light line (dashed). The result 
agrees well with the finite array dispersion (blue open circles 
taken from (A)), but not with the quasistatic approximation 
(thin curves, dotted for longitudinal mode). (C) Imaginary 
part of the frequency for the modes in (B). The damping rate 
diverges near the light line for transverse modes. Inset in (C): 
Propagation lengths (amplitude 1 /e lengths) versus frequency 
for transverse (red solid circles) and longitudinal modes (green 
open squares). The horizontal dashed line indicates the inter- 
particle spacing d. 



